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Intersections  et  Convergence 


Paul  de  Faget  de  Casteljau 


Abstract.  C’est  en  1958,  que  j’ai  rejoint  l’equipe  de  fraisage  numerique, 
CAO  de  la  prehistoire,  animee  par  Monsieur  de  la  Boixiere  : Farfelus 
pour  les  uns,  fous  pour  d’autres,  personne  en  dehors  de  notre  chef  n’aurait 
ose  parier  un  seul  kopek  sur  l’avenir  de  cette  technique.  Les  problemes 
a surmonter  etaient  aussi  nombreux  que  fort  ardus  : Par  exemple,  celui 
d’intersections,  ou  racines.  Nos  ordinateurs  se  bloquaient.  De  toute  ur- 
gence,  il  fallait  mettre  sur  pied  une  methode  efficace,  en  harmonie  avec  les 
exigences  de  ce  nouvel  outil.  Apres  tant  d’essais  d’autant  plus  volumineux, 
que  longtemps  infructueux,  comme  tout  parait  simple,  quand  on  presente 
le  joyau  final!  Mieux  encore,  une  idee  directrice,  qui  comblera  d’aise  les 
ferus  de  fractions  continues,  nous  y conduit  en  toute  logique. 


§1.  Un  Peu  d’Histoire 

Pour  demarrer  l’etude  des  formes  de  C.A.O.,  quelques  instants  de  reflexion 
pouvaient  sufflre,  pour  se  convaincre  que,  pour  concilier  les  imperatifs  du 
dessin  a ceux  d’une  production  de  calculs,  la  seule  solution  acceptable  restait 
l’utilisation  en  geometrie  affine  des  formes  parametriques  polynomiales  a tra- 
vers  l’algorithme  devenu  celebre.  La  necessity  d’obtenir  une  solution,  unique 
pour  tout  parametre,  interdisait  des  calculs  aventureux.  Meme  si  le  calcul  de 
fonctions  du  genre  arc  cos(®),  au  voisinage  de  x = 0,  aurait  gagne  du  point  de 
vue  logique,  de  l’utilisation  de  formes  y/x,  deja  multiforme,  les  mathematiques 
en  sont  restees  aux  series  entieres!  (rarement,  les  fractions  rationnelles,  frac- 
tions continues...). 

II  n’en  restait  pas  moins  a traiter  de  nombreux  problemes,  qui,  chacun, 
pouvait  me  donner  l’impression  d’avoir  a passer  le  concours  de  l’agregation, 
me  creant  des  cauchemars  nocturnes,  tout  en  regrettant  amerement  le  bon  et 
solide  poste  de  professeur  qui  aurait  pu  etre  le  mien. 

“M’sieur,  ga  ne  marche  pas!”,  phrase  rituelle,  qui  m’annongait  encore  une 
nouvelle  catastrophe,  a cette  epoque  heroi'que,  pour  ne  pas  dire  prehistorique. 
Cela  m’arrachait  a mes  meditations  et  reveries  ou  je  me  figurais  dans  une 
sorte  de  paradis  terrestre,  ou  les  memes  visiteurs  me  felicitaient  de  ma  derniere 
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trouvaille  d’un  “Ce  n’est  vraiment  pas  mal,  votre  true!” , phrase  toujours  restee 
dans  le  virtuel  de  ma  pensee. 

Parmi  ces  innombrables  problemes,  qui  me  ramenaient  brutalement  au 
sens  des  realites,  il  y eut  le  probleme  des  intersections.  La  encore,  pas  ques- 
tion d’adopter  une  attitude  universitaire,  autour  d’une  savante  discussion  sur 
l’existence  et  la  nature  d’eventualites  favorables.  II  fallait  “produire”  une 
racine  et  une  seule,  et  de  surcroit,  de  valeur  acceptable.  Cette  fois,  l’unicite 
devenait  difficile  a garantir,  surtout  au  voisinage  de  racines  multiples,  voisines 
ou  confondues.  Aucune  des  tentatives  de  tenir  compte  des  termes  non  lineaires 
(Whittaker,  etc...)  pour  ameliorer  la  methode  de  Newton  ne  s’est  vraiment 
imposee.  Et  les  exploitations  en  ordinateurs  atteignaient  la  limite  du  quart 
d’heure  qui,  par  precaution,  en  limitait  la  duree!  Et  de  plus,  il  fallait  agir 
vite,  tres  vite...pour  reprendre  un  probleme  tres  classique  a zero,  en  poussant 
notre  “cri  de  guerre”  : “Mais,  qu’alors  y faire?”  (jeu  de  mots:  allusion  a un 
calorifere). 


§2.  Expose  du  Probleme 

On  veut  calculer  la  valeur  de  x,  qui  verifie  l’equation 


P(x)  = 0, 


oil  P(x)  est  un  polynome  (ou  un  developpement  limite,  en  serie  de  Taylor, 
d’une  fonction 

Ce  probleme  se  pose  aussi  a plusieurs  inconnues;  ainsi,  pour  n = 3,  calcul 
de  x,  y,  z pour  que 

( P(x,y,z)  =0, 
lQ(x,y,z)  =0, 

{ R(x,  y,  z)  =0. 

Pour  le  cas  n = 2,  j’ai  presente  au  congres  de  SchloC  Dagstuhl  (Saarland 
1996)  une  fort  jolie  methode,  habile  manipulation  de  geometrie  projective,  tres 
efficace;  malheureusement,  tous  mes  efforts  de  l’etendre  a d’autres  valeurs  de 
n ont  echoue.  J’ajouterai  a l’adresse  des  chercheurs,  qui  auraient  la  temerite 
d’etudier  le  cas  n — 3,  que  j’avais  bien  cerne  la  difficulte  qui  rendait  ce  cas 
impossible. 

Les  seules  equations  qui  conduisent  a une  solution  unique,  sont  du  premier 
degre.  Il  nous  reste  la  possibilite  de  tenter  une  identification  autour  de  x = 0 
du  polynome  P{x)  a une  fraction  rationnelle  dont  le  numerateur  soit  lineaire. 


n = 1 
n = 3 


P{x)  ~ 


P(x,y, 


a + bx 

1 + dx  + ex2  + ... 

. a + bx  + cy  + dz 
1 + 5^  hpxayPz~< 


On  en  deduit  une  suite  d’approximations  de  la  racine,  chacune  utilisant 
les  precedentes,  de  faqon  optimale.  Cela  peut  se  prouver  par  l’analyse,  et 
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le  verifier  sur  des  exemples  concrets : elle  est  la  meilleure,  comparee  a toute 
manipulation  du  meme  type. 

L’identification  brutale,  terme  a terme,  conduit  a 

A + Bx  + Cx2  ~ 

A + Bx  + Cx 2 + Dx3  ~ 

A + Bx  + Cx 2 + Dx 3 4-  ExA  ~ 


AB  + ( B2-AC)x 
B — Cx  + . . . 

A(B2-AC)  + [A(AD-BC)  + B{B2-AC)]x 
B2~AC  + {AD-BC)  x + ( C2-BD)x 2 + . . . 
num 


num  = A[A(AD  - BC)  + B(B2  - AC)] 

+ [A2(BD  - AE)  + AB(AD  - BC)  + ( B 2 - AC)2]x 
den  = [A(AD  - BC)  + B{B2  - AC)\ 

+ [A(BD  - AE)  - C(B2  - AC)\x 

+ [B(AE  - BD)  - C(AD  - BC)\x2 

- [A{D2  - CE)  + B{BF  - CD)  + C(C 2 - BD)] x3  + . . . 

Seuls  les  numerateurs  nous  interessent: 


B 

A 

0 

0 

0 

C 

B 

A 

0 

0 

D 

C 

B 

A 

0 

E 

D 

C 

B 

A 

F E D C B 


Si  on  appelle  Ui  la  succession  des  determinants  utilisant  les  i premieres  lignes 
et  colonnes  de  la  matrice  ci-dessus,  en  posant  Uq  = 1,  on  obtient  la  sequence 
d’approximations 

AUi+ 1 

** = uT' 

x\  — —jj  (Newton),  x%  coincide  avec  Whittaker,  mais  non  X3  et  au  dela. 
Mais  il  y a mieux  a faire  ! 


§3.  Algorithme 

Une  etude  complete  de  cette  question  prouve  que  l’on  obtient  le  meme  resultat 
en  utilisant  un  algorithme  qui  rappelle  le  schema  de  Horner. 

II  s’agit  de  calculer  la  racine  x,  en  partant  de  la  valeur  x = 0 d’une  serie 
entiere  type  MacLaurin  ecrite  sous  la  forme 


a0  = ai_x  + a2x2  + . . . + anxn  + ... 


On  ecrit  alors  l’algorithme  sous  la  forme : 


12 


P.  de  Casteljau 


clq  = aiN 

donne  N ~ Newton 

a-o  = (a  i + a2N)P 

” P ~ Whittaker 

1 

no  = [ai  + (a2  + a3  N)P]Q 

” Q ~ Fractions  continues 

22 

ao  = {ai  + [°2  + (ffl3  + n4N)P]Q}/i 

” R 

53 

130' 

. . . etc  . . . 

” S,T,... 

On  a indique  a droite  le  nombre  de  possibilites  d’utilisation  des  approxi- 
mations anterieures,  dans  un  ordre  different,  ou  avec  repetition.  Mes  etudes, 
a l’epoque,  m’ont  conduit  a prouver  qu’aucune  ne  fait  mieux,  avec  les  hy- 
potheses adoptees  (le  nombre  130  est  la  somme  des  trois  precedents,  pour 
tenir  compte  de  toutes  les  approximations  nouvelles). 

Cette  methode  est  d’autant  plus  efficace  que  la  sene  a.  inverser  est  plus 
convergente  ; cela  se  verifie  pour  e1,  sin(x),  sh(x),  pour  calculer  ln(x),  Arc 
sin(x),  Arg  sh(x). 

On  va,  a titre  d’exemple,  calculer  le  logarithme  neperien  de  2,  la  ou  la 
serie  ln(l  + x),  au  bord  de  son  cercle  de  convergence,  est  particulierement 
desesperante. 

II  convient  de  n’y  voir  qu’un  exemple,  et  non  une  invitation  a faire  le  calcul 
des  logarithmes  de  cette  faqon-la.  Notons  encore  tout  l’interet  d’inverser  les 
lignes  trigonometriques  par  l’utilisation  simultanee  du  sinus  et  du  cosinus,  ou 
encore  pour  un  triangle  rectangle  les  cotes  S,C,H  du  sinus,  du  cosinus  et 
l’hypotenuse. 

Pour  les  lecteurs  interesses  par  cette  question,  on  peut  demontrer  deux 
sequences  de  reduites  de  fractions  continues,  co'incidant  une  fois  sur  deux,  ou 
encadrant  par  exces  et  par  defaut,  la  valeur  exacte  de  la  fonction  : l’une  de 
ces  sequences  se  rattache  au  developpement  de  Arc  tg(ip/2)  = t avec  S — 
2 1,  C = 1 — f2,  H = 1 + t2,  l’autre  est  originale.  A l’epoque  ou  le  calcul 
des  fonctions  n’etait  pas  encore  tres  au  point,  au  niveau  des  ordinateurs,  les 
formules  simplifiees  de  Arc  sin  et  Arc  cos  nous  ont  rendu  un  grand  service. 
Ainsi,  tp  # S.  qui  n’est  pas  ridicule  a 5 = 1,  C = 0 : Mais 

cela  n’est  pas  le  sujet  de  cet  expose. 

Exemple  numerique : 

Appliquons  cet  algorithme  au  calcul  de  x,  defini  par : 

2 = e-  ~ 1+  |(1  + |(1  + |(1  + . . . . (1  + |)  ...  )...))) 


On  trouve  ainsi  : 


on  pourrait  encore  considerer  les  possibilites  obtenues,  en  prenant  les 
formes  polaires  des  N,P,Q.  Ainsi,  a la  seconde  ligne,  au  lieu  de  a\P,  on 
ecrirait  ai  — , etc... 
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1.00 
0.666  6 
0.692  30 
0.693  333 
0.693  1608 
0.693  145  41 
0.693  146  977 
0.693  147  196  49 
0.693  147  183  337 
0.693  147  180  436 


0.693  147  180  559  945  309 

soit  environ  un  gain  d’une  decimale  a chaque  iteration,  ce  qui  reste  plus 
qu’honorable  compare  avec  les  series  traditionnelles  ln(l  + u)  avec  u = 1 
ou  ln(ji^)  avec  u — |,  bien  plus  convergente. 

§4.  Promenade  au  Pays  du  Logarithme 

II  est  toujours  passionnant  de  reprendre  a son  compte  l’etude  des  fonctions 
elementaires,  qui  conduisent  a des  developpements  d’une  telle  richesse,  qu’ils 
paraissent  inepuisables.  Ainsi  en  est-il  avec  le  logarithme,  l’etude  qui  suit, 
et  qui  n’est  qu’un  tout  petit  apergu  des  proprietes  que  revele  l’analyse,  sans 
parler  des  proprietes  devoilees  par  l’arithmetique. 

L’etude  precedente,  reprise  de  fagon  litterale,  founit  la  succession  d’appro- 
ximations 

(x  _ - 2 • 3 (^-!)(^+1)  . a (x-l)(x2+4s+l)  . 5(x-1)(x+1)(x2+10x+1)  . 

1 ) > ^(x+l)  ’ ° x2+4x+l  ’ *(x+l)(x2  + 10x+l)  > ° x4+26x3+66x2+26x+l  ’ 

On  constate  deja  plusieurs  regies  empiriques.  Sans  insister,  posons  : 
si  n = 2p  : sin  = 2p+l:  n (x-i)(x+i)%|^. 

On  peut  encore  reduire  l’ecriture  de  ces  polynomes  reconnus  symetriques,  en 
posant  x + A — 2 = 2t  et  si  kq  = 1 pour  les  polynomes  P et  kq  = q pour  les 
polynomes  Q : 

Q ou  P = 1 

Pi  ou  Qi  = t + 3 

P2  ou  Q2  =t2  + 15 k2t  + 30&3 

P3  ou  Q3  — t3  + 6Zk2t2  + 420&3<  + 630A:4 

P4  ou  Qi  = t4  + 255 k2t3  + 4410 k3t2  + 18900fc4<  + 22680fc5 


On  peut  encore  comparer  ces  fonctions  avec  les  reduites  du  developpement 
de  ln(a;)  en  fractions  continues,  tres  inhabit uelles,  parce  qu’elles  sont  centrees 
autour  de  x = 1 au  lieu  du  zero  habituel.  Elies  sont  tres  faciles  a obtenir, 
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par  integration  de  formules  Q(a;)ln(a:)  — P(x)  = k(x  — l)n,  alternees  avec  des 
changements  de  x en  l/x,  qui  change  le  signe  de  ln(a:). 

On  remarquera  les  carres  des  coefficients  du  binome,  dans  les  expressions 
des  denominateurs  de  rang  pair  (cela  se  demontre). 

. 2(x-l)  . x2+4x-5  . 3(x2-l)  . x3  + 18x2-9x-10  . 

> x+1  ’ 4x+2  1 x2+4x+l  ’ 9x2  + 18x+3  ’ 

11x3+27x2-27x-11  . 3x4  + 128x3  + 108x2-192x-47  . 25x4  + 160x3-160x-25 

3(x3+9x2+9x+l)  > 12(4x3  + 18x2  + 12x+1)  > 6(x4  + 16x3+36x2  + 16x  + l) 

Si  l’on  pose  Dn  = en  utilisant  la  valeur  des  reduites  paires, 

on  peut  meme  arriver  a etablir  cette  formule  particuliement  brillante  : 


i \ or^>  , diD2  , D2D3  D3Di 

ln( *) ) ~ 2 [ T + — + — + — + • • 


■ + 


D-n-lDn 


Remarque.  L’etude  arithmetique,  a laquelle  j’ai  fait  allusion,  obeit  a une 
demarche  opposee  a celle  de  la  table  des  nombres  premiers.  On  ne  retient 
des  nombres  entiers  que  ceux  dont  les  facteurs  premiers  ne  sont  que  les  n 
premiers  nombres  de  cette  table  des  nombres  premiers.  Les  inegalites  sur  leurs 
logarithmes  qui  en  decoulent  fournissent  une  serie  de  solutions  aux  problemes 
de  Dirichlet-Hermite.  On  trouve  ainsi  une  explication  rationnelle  aux  gammes 
musicales...  Voir  a ce  sujet  les  derniers  chapitres  de  mon  livre  “Quaternions”. 

Ces  considerations  s’appliquent  aussi  a des  fonctions.  Ainsi  trouver  les 
meilleures  solutions  polynomiales  a un  degre  donne,  sensiblement.  proportion- 
nelle  a un  arc,  et  a leurs  sinus  et  cosinus.  On  peut  y rajouter  la  tangente, 
etc... 


§4.  Cas  de  Plusieurs  Variables 
Le  regard  de  l’analyse 

Le  probleme  est  trop  vaste  pour  etre  expose  en  quelques  lignes.  Nous 
renvoyons  aux  ouvrages  specialises,  pour  constater  qu’il  ne  suffit  pas  de  se 
donner  p points,  pour  definir  la  valeur  de  p coefficients  d’une  interpolation  a 
plusieurs  parametres.  II  intervient  alors  la  notion  de  silhouette  du  polynome 
associe  a une  grille  de  points  donnes.  Les  silhouettes  les  plus  considerees,  a 
deux  variables,  peuvent  etre  ainsi  triangulaires  J2  x°y0  avec  a + /3  < n,  ou 
encore  rectangulaires  si  a < n,  (3  < q et  carrees  si  n = q.  Paradoxalement, 
c’est  avec  les  grilles  les  plus  regulieres  qu’on  rencontre  le  plus  de  contraintes 
a respecter. 

C’est  en  fonction  de  la  notion  de  sihouette  que  Ton  imposera  un  certain 
ordre  de  classement  aux  termes  d’un  polynome  a plusieurs  variables.  Pour 
determiner  une  fraction  rationnelle  a numerateur  lineaire  au  voisinage  d’un 
point,  on  peut  imaginer  l’ecriture  d’un  determinant  dont  la  premiere  ligne 
s’ecrit : 

i,  xi  y,  f,  xf,  yf > x2f,  xyf, 

la  deuxieme  ligne  exprime  que  / passe  par  le  point  X,  Y,  F (determinant  nul,  si 
deux  lignes  egales)  ensuite,  un  point  voisin  en  X (puis  en  Y)  et  par  difference 
avec  cette  seconde  ligne  s’exprimera  a partir  des  derivees  part.ielles  en  X (puis 
en  Y).  On  poursuivra  les  lignes  suivantes  avec  les  derivees  d’ordre  superieur. 
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Mais  voila  ! On  se  heurte  a un  mur:  le  probleme  est  indetermine. 

“Qu’alors  y faire”  ? Faut-il  renoncer  ? Le  probleme  exige  une  solution, 
brillante,  ou  batarde,  peu  importe,  a condition  d’en  trouver  UNE  ! Demarche 
industrielle,  et  non  universitaire. 

Alors,  miracle  ! Quelle  que  soit  la  faqon  dont  on  leve  l’indetermination, 
par  un  choix  arbitraire  d’un  terme  de  chaque  degre  du  denominateur,  on 
retrouve  le  meme  numerateur;  or  c’est  de  lui  et  de  lui  seul  qu’on  a besoin,  et 
tout  rentre  dans  l’ordre. 

Ici,  honnetement,  je  dois  avouer,  que  je  n’ai  developpe  ce  calcul,  assez 
fastidieux,  que  dans  le  seul  cas  de  deux  variables. 

Retour  a l’algorithme 

II  est  plus  rapide  de  faire  un  raisonnement,  que  d’aucuns  jugeront  fort 
acrobatique,  mais  parfaitement  correct.  Dans  chaque  equation,  on  suppose 
le  probleme  resolu,  pour  toutes  les  variables,  excepte  une,  d’abord  x (puis 
y,z,...).  On  applique  alors  notre  algorithme  au  calcul  de  la  premiere  ap- 
proximation de  Newton  done  Nx  (puis  Ny ,NZ).  On  resoud  le  systeme,  ce  qui 
donne  Nx,  Ny,  Nz. 

On  utilise  ces  valeurs  approchees,  comme  on  sait  le  faire  pour  determiner 
Px  (puis  Py,Pz)  et  par  suite  Px,Py,Pz.  Une  troisieme  etape  permettra  le 
calcul  de  Qx,  Qy,  Qz,  etc  ... 

La  meilleure  silhouette  qui  convienne  est  done  carree  {ot  < n,  (3  < n,  7 < n), 
ou  completee  par  des  0. 

Remarque:  Quel  que  soit  le  nombre  de  variables,  l’optimisation  du  calcul 
montre  que  l’ideal  se  situe  entre  la  seconde  iteration  Px , Py , Pz , et  la  troisieme 
Qx,  Qy,  Qz-  Apres  il  est  preferable  de  revenir  a un  developpement  de  Tay- 
lor centre  sur  ce  nouveau  point,  sauf  cas  de  singularite  extraordinaire,  qu’en 
principe  on  ne  doit  pas  rencontrer  dans  la  pratique  industrielle. 

Voila,  en  raccourci,  l’expose  de  ce  probleme  tres  utile.  Son  etude  ap- 
profondie  pourrait  faire  l’objet  d’une  these  complete.  La  discussion  semble 
ici  escamotee:  ce  n’est  pas  tout  a fait  vrai;  car  si  une  discussion  s’impose, 
c’est  que,  dans  la  pratique  industrielle,  il  n’y  a plus  de  solution  pratique.  Par 
exemple,  on  ne  cherche  jamais  l’intersection  d’une  forme  avec  une  verticale, 
qui  coupe  la  forme  en  trois  points,  interieurs  au  carreau  considere.  Comment 
l’usinerait-on  ? Comment  l’emboutir  ? 
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